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ABSTRACT 


The  dynamics  of  the  coupling  of  linear  internal  gravity 
< waves  and  linear  surface  gravity  waves  on  the  ocean  is  studied 

using  a Hamiltonian  formalism  and  action-angle  variables.  The 

. V 

dynamic  equations  are  solved  both  numerically  and  in  some  anal- 
ytic approximations.  The  results  compare  favorably  with  the 
interaction  experiments  of  Lewis,  Lake  and  Ko  the  "resonant 
triad"  experiments  of  Joyce  and  some  satellite  observations 
of  Apel  et  al.  The  growth  time  for  internal  waves  generated 
by  the  resonant  interaction  of  surface  waves  is  calculated 
using  the  Garrett-Munk  ocean  model  and  the  Phillips  spectrum 
for  surface  waves.  Energy  exchange  rates  are  deduced. 
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I.  INTRODUCTION 

The  possible  importance  of  surface  waves  for  the  genera- 
tion of  internal  waves  is  mentioned  by  Phillips  (1966)  and 
models  for  this  have  been  studied  by  Ball  (1964) , Thorpe  (1966), 
and  Kenyon  (1968) . Tank  experiments  to  study  the  interaction 
of  surface  and  internal  waves  have  been  reported  by  Lewis, 

Lake,  and  Ko  (1974)  and  by  Joyce  (1974) . 

Synthesis  of  observations  of  internal  waves  in  the  deep 
ocean  have  been  made  by  Garrett  and  Munk  (1972a) , (1975) , who 
have  proposed  an  explicit  equilibrium  spectrum  for  internal 
waves.  Suggested  sources  of  energy  to  drive  internal  waves 
are  tidal  currents,  atmospheric  pressure  and  stress  fluctuations 
and  surface  waves  £see  for  example  Phillips  (1966) , Chapter  5 
or  Thorpe  (1975) j . The  relative  importance  of  these  mechanisms 
at  various  internal  wavenumbers  is  not  presently  known.  If  it 
is  not  saturated  one  might  anticipate  that  energy  is  fed  into 
the  internal  wave  spectrum  at  longer  and  intermediate  wave- 
lengths and  cascades  by  nonlinear  interactions  to  the  shortest 
internal  wave  wavelengths,  where  it  is  dissipated  by  wave 
breaking. 


The  purpose  of  this  paper  is  to  study  the  generation 
of  internal  waves  by  surface  waves.  We  shall  simplify  the 
dynamics  by  treating  the  surface  wave  and  internal  wave  motion 
in  the  linear  approximation.  The  nonlinear  coupling  will  be 
treated  in  the  lowest  non-vanishing  order.  To  further  simplify 
our  analysis  we  shall  assume  (not  unarbitrarily ) that  the 
wavenumbers  and  frequencies  of  the  surface  waves  studied  are 
large  compared  to  the  wavenumbers  and  frequencies  of  the 
internal  waves  of  interest.  The  resulting  dynamics  will  be 
expressed  in  Hamiltonian  form  and  the  Hamiltonian  equations 
of  motion  expressed  using  action-angle  variables.  Both  analytic 
approximations  and  numerical  integration  of  these  equations 
will  be  described. 

For  the  purposes  of  this  paper  we  shall  model  the  ocean 
as  follows.  In  the  absence  of  wave  motion  the  surface  is 
assumed  to  locally  coincide  with  the  plane  z = 0 of  a rectan- 
gular coordinate  system.  The  bottom  is  assumed  to  coincide 
with  the  plane  z = -H.  In  order  to  use  discrete  Fourier 
expansions,  a large  rectangular  ocean  area  Aq,  with  periodic 
boundary  conditions,  will  be  considered. 

The  density  p(z)  is  assumed  to  have  the  following 
characteristics : 

p(z)  = p , a constant,  for  - D < z < 0. 


A thin  thermocline  with  density  change  6p  occurs  at  z = - D. 

For  -H  < z < -D,  p(z)  is  a monotonically  decreasing  function  of 
z . 

It  is  the  Brunt-Vaisala  frequency  N(z)  which  will  appear  in 
our  equations  for  the  internal  waves  £see,  for  example,  Phillips 
(1966),  Chapter  5]  . The  chosen  model  for  the  ocean  implies 
that 


N ( z ) =0  , -D  < z < 0, 

N (z)  = £~(g/P0)  , -H  < z < -D 


and 


-D+6 

N2 (z)dz  = g 6p/ PQ  , (1) 

where  (26)  is  a small  interval  somewhat  greater  than  the 
thickness  of  the  thermocline.  The  quantity  g represents 
the  acceleration  of  gravity. 

Although  our  dynamical  equations  will  be  obtained  for 
this  rather  general  model,  the  numerical  calculations  will  be 
based  on  either  of  two  more  specific  models.  The  first  is 
that  of  Garrett  and  Munk  (1972a) : 

N (z)  =0  , -D  < z < 0, 

N(z)  = Nq  exp  £(z  H D)/bj,  -H  < z < -D, 
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N = 5.2x10  ^ sec  \ 
o 

H = 5000  meters 

b = 1200  meters.  (2) 

The  internal  wave  eigenmodes  for  this  model  will  be  obtained 
and  used  in  the  WKB  approximation. 

The  second  model  which  we  shall  use  is  the  somewhat 
unrealistic,  but  convenient,  thin  thermocline  model: 

N(z)  = 0,  except  near  z = -D , 

-D+S 

N2(z)dz  = g Sp/pQ  . (3) 

A constant  Brunt-Va'isalci*  frequency  model  is  sometimes  used. 
Kenyon  (1968) , for  example,  has  studied  the  coupling  of  sur- 
face and  internal  waves  using  this  model.  We  shall  see  that 
this  model  leads  to  very  much  weaker  coupling  than  the  Garrett- 
Munk  model  or  the  thin  thermocline  model  ^already  observed  by 
Joyce  (1974) J and  is  presumably  unrealistic  for  studying  the 
generation  of  internal  waves  by  surface  waves. 

The  derivation  of  the  equations  for  coupled  surface 
and  internal  waves  will  be  outlined  in  Section  II.  Applica- 
tions to  the  tank  experiments  of  Lewis  et  al.  (1974)  and  Joyce 
(1974)  will  be  described  in  Section  III.  Miscellaneous  appli- 
cations, including  the  generation  of  internal  waves  by  a large 
amplitude  swell  [Apel  et  al.  (1975a)] , will  be  described  in 
Section  IV.  The  rate  of  generating  internal  waves  by  a 
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saturated  spectrum  of  surface  waves  will  be  studied,  including 
the  dependence  on  mixed  layer  thickness  D and  thermocline 
strength  6p/pQ  in  Section  V.  The  results  of  some  numerical 
calculations  of  the  development  of  a mottled  appearance  of 
the  surface  as  the  internal  waves  grow  in  amplitude  ^Apel,  et 
al.  (1975b)]  will  be  shown  in  Section  VI. 


II.  THE  coupling  of  surface  and  internal  waves 

The  equations  describing  mutual  coupling  of  surface 
and  internal  waves  will  be  presented  in  this  section.  In  the 
mixed  layer,  corresponding  to  -D  < z < 0,  the  flow  is  treated 
as  irrotational  and  the  fluid  velocity  u can  be  expressed  as  the 
gradient  of  a potential  <p(r,z).  ^Here,  and  elsewhere  in  this 
paper  we  let  r = (x,y)  represent  a vector  in  the  horizontal 
plane  z = 0.] 

We  have  assumed  that  the  internal  wave  modes  of  inter- 
est to  us  have  much  lower  characteristic  frequencies  and  much 
longer  wavelengths  than  have  the  surface  waves.  This  suggests 

A 

writing  the  velocity  potential  <(>  in  the  form 

<f>(r,z)  = <j>  (r , z)  + $ (r,z)  . (4) 

Here  <J>  contains  the  high  frequency,  high  wavenumber  components 

A 

(surface  waves)  of  <j>  and  contains  the  low  frequency,  low 

A 

wavenumber  components  (internal  waves)  of  <j> . Similarly,  we 
can  write  the  vertical  displacement  h,  due  to  wave  motion,  of 
the  ocean  surface  as  a sum  of  a high  frequency  part  £ and  a 
low  frequency  part  H: 

h(r,t)  = ;(r,t)  + H (r , t)  . (5) 

Bernoulli's  equation  and  the  kinematic  boundary  condition 
on  the  surface  z = h(r,t)  are,  respectively, 
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z = h 


r 


II  + I (vj)2  + gh  = 0 , 

A 

V£  + (V>-<V>  ’ If  ' 2 = h ' (6) 

where  V is  the  two-dimensional  gradient  operator,  i.e., 
s 

vs  ~ \3x  ' 3y } . 

For  small  amplitude  waves  the  fluid  motion  below  the 
ocean  surface  is  described  by  the  equation  J see  for  example, 
Phillips  (1966),  Chapter  5J 

■ilf-L.  + V2\  w(r,z,t)  + N2(z)  V2  w(r,2,t)  =0.  (7) 

3t2  \3z2  s)  ' s ~ 

Here  w is  the  vertical  component  of  the  fluid  velocity  and  N 
is  the  Brunt-VSTisalk*  frequency  given  by  Eq.  (1)  . We  have 

3^, 

assumed  that  N = 0 in  the  mixed  layer,  where  w = g--  . In  the 
region  of  the  thermocline  and  below,  we  assume  that  the  flow 
w is  associated  with  internal  waves  only. 

Because  we  have  supposed  the  wavelengths  of  the  surface 
waves  to  be  small  compared  to  those  of  the  internal  waves,  and 
also  not  significantly  larger  than  the  thickness  D of  the 
mixed  layer,  the  nonlinear  terms  responsible  for  the  coupling 
between  the  surface  and  internal  waves  are  contained  in  the 
surface  boundary  conditions  (6) . 


Near  the  surface,  the  low  frequency  fluid  velocity 
associated  with  internal  waves  is 

U(r,t)  =V$|z=0  = U(r,t)  + k U2(r,t)  (8) 

where  U and  Uz  are  the  horizontal  and  vertical  currents, 
respectively.  The  corresponding  high  frequency  velocity  field 
associated  with  surface  waves  is 


The  properties  of  internal  waves  permit  us  to  assume  that 

|uz|  « |y| 

N « |c| 

|V*(H)/H|  « |?^(5)/S|  . (10) 

The  first  two  inequalities  describe  the  fact  that  for  internal 
wave  motion  there  is  little  displacement  of  the  ocean  surface. 
The  third  inequality  states  that  the  surface  curvature  due  to 
internal  waves  is  expected  to  be  much  less  than  that  due  to 
surface  waves. 

To  describe  the  surface  wave  motion,  we  first  define 
the  high  frequency  velocity  potential  on  the  surface  z = H(r,t) 
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(11) 


<J>  _(r,t)  = 4>  ( r , z = H(r,t),  t) 

and  extract  the  high-frequency  pa). t of  Eqs.  (6)  . With  use  of 
the  inequalities  (10) , we  easily  obtain  to  the  specified  order 
of  non-linearity,  the  coupled  equations 


(&  + 2‘7s)*s  + « = o , 

{jt  + 2‘7s)  « + « V"  - ®»s  <12> 

Here  we  have  introduced  the  operator 

® = (-V‘)  (13) 

and  understand  that  it  will  always  act  on  a Fourier  series. 
Laplace's  equation,  V 4>=0,  has  been  used  to  write 

H-  (7sH»-7s*]z=H  +<9r[(7sH)*] 

in  the  second  of  Eqs.  (12) . 

We  next  extract  the  low  frequency,  long  wavelength 
portion  of  Eq.  (6) . On  following  an  argument  of  Phillips, 
[l966,  Eq.  (5.2.12)j  we  obtain  with  a little  straightforward 
algebra  the  surface  boundary  condition 


Here  we  have  written  <^r)>L  as  representing  the  "low  frequency 
long  wavelength  part"  of  the  quantity  r. 

To  continue,  we  suppose  that  in  the  absence  of  surface 
waves  a complete  set  of  eigenmodes  for  Eq.  (7)  are  known  to 
be  of  the  form 

w(r,z,t)  = W°;L  (z,  exp£[t.r  - n(j,L)t)J  , ,15) 


where 


Wi  T (0)  = w°  r (-«)  = 0 
J » ■“  J , L 


Thus,  (j,L)  describes  a free  internal  wave  of  mode  number 
j Jj  = 1,2, ...in  the  notation  of  Garrett  and  Munk  (1972a)] 
horizontal  wave  number  L and  frequency  fi(j,L).  The  normaliza- 


tion of  W.  is  so  chosen  that 
D > ij 


[W°,L  <2»] 


= L' 


z = 0 


(16) 


Here,  and  in  the  following,  we  write 


F',z>  s H - 


etc . 

Equations  (7)  and  (16)  imply  that  in  the  mixed  layer, 
where  N=0,  we  can  write 
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w.  T (z)  = - l sinh  (Lz) 
J 


t 


-D  < z < 0 . 


(17) 


At  the  thermocline  W _ 

D /L 


is  continuous,  but 


above 


- (MU 


below 


2,,o 


wj  fL(-D)/n‘(j,L) 


Finally,  in  the  range  -H  < z < -D,  W?  L is  found  by  integration 
of  Eq.  (7) . 

The  above  expression  at  the  thermocline  results  in  the 
dispersion  relation 


(j  = [(g  <5p/po)  ,L)  + coth(LD)] 


(18) 


where 


(J\  (j,L)  = L-1(in  W° /L(z)) 


evaluated  for  z just  below  the  thermocline. 

To  obtain  the  corresponding  solutions  to  Eq.  (7)  with 
the  surface  boundary  condition  (14) , we  write  (in  the  rectan- 


gular area  A of  ocean  surface  being  considered) 


^(r,z)  = J]  4>  (z)  exp(i  L*r)  , -D  < z < 0 , (19) 

j , L J'" 


and 


w(r,z)  = J W.^L  (z)  A.  T (t)  exp(i  L-r) 


j > 1* 


j/L 
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(20) 


A 


■MP'J  ' J.  ■■  .."•'IIWMP4.  IHUIH  ^yinp  mmw*  W'l  w > 


The  form  of  Eq.  (20)  permits  us  to  so  normalize  W.  T that 

J 


(21) 


On  writing,  from  Eq.  (14) , 


V <«7s  *S>L  £ £ rL  exp(i  ’ 

and  on  substituting  this  and  the  expansion  (20)  into  Eq.  (14) , 
we  obtain 


(0)  A. 


L 


(22) 


Again,  if  we  substitute  the  expansion  (20)  into  Eq.  (7) , and 
(for  the  moment)  drop  the  (j,L)  labels,  we  obtain  the  equation 


w"  + L2  [n2/A2  - lj  W = 0 


(23) 


A 


2 


f 


where  A E ^ 


etc. 


We  now  consider  to  be  small,  ie  .,  the  vertical 
velocity  of  the  internal  wave  is  considered  small  with  respect 
to  the  vertical  variation  of  the  free  surface,  and  write, 

(24) 
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w = w~  + 6w  , 


A 2 = a2  + 6A2  ' 


where  W°  represents  one  c ” the  eigenmodes  (15)  and  is  one 
of  the  eigenfrequencies  (18).  Equation  (22)  can  be  re-expressed 
in  the  form 


6WA  = r 


(ie  MWJ=0atz=0), 

To  obtain  an  equation  for  A(t) , we  consider  the  integral 


o 

J = J w[w"  + L2(N2/A2  - l)w]  dz  = 0 


This  integral  is  varied  using  F.q.  (24)  and  can  be  written  as 
j = J°  + 6j,  where  J°  and  6J  are,  respectively,  of  order  zero 
and  one  in  the  small  quantity  T.  After  some  manipulation,  the 
equation  6j  = 0 can  be  put  into  the  form 

° 2 1 _1 

a + o2a  = |n4r(w°(z)VL  «1  [l2  f n2(z>  (w°(z))  dz  » (27) 


[a4r(w°<z))'l2J 

|l2  r n2(Z) 

(w°  ( z )^  dzj 

J 

1 J 

l -H 

valid  to  first  order  in  Y. 


To  simplify  the  right-hand  side  of  Eq.  (27) £and  putting 
back  the  (j,L)  labelsj it  is  convenient  to  define  the  quantity 
g'(j,L)  with  the  relation 


(28) 


J N ( z ) L(z)J  dz 
-H 


g' ( j ,L) 


= L2  g' ( j ,L)  sinh2  (LD) 

Equation  (17)  was  used  in  the  final  step  here.  The  expressions 
(16),  (21)  and  (28)  allow  us  to  put  Eq.  (27)  into  the  form 


* , t (0)  + n 2 ( j , L)  <D.  . (0) 

3 / 1 D • " 

= fi4(j,L)  rL  /[L2g  '(j,L)  sinh2  (LD)J 

= C ( j , L)  Tl  . (29) 

We  can  formally  return  to  the  r coordinate,  defining 
[see  Eq.  (19)j 

<t>.(r)  = E (o)  exp(iL*r) 

3 ' L j,L 


Then  Eq.  (29)  can  be  re-written  as 

*.(r)  + fi2(j,d»  *j(r)  = C(j,®)Vg-  (30) 

] L 

Here  ft  and  C are  the  appropriate  operators  to  give  Eq.  (29) 

when  the  modal  expansion  is  introduced. 

The  vertical  displacement  -(r,t)  at  the  thermocline  is, 
in  the  linear  approximation , obtained  on.  integrating  the 
equation 


" -=  Z 

j 


E A-t  T W?  t exP(i  L-r) 

j , l 3't?  J 'h  - - 
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= - T L sinh  (LD)  $.  T (0)  exp(i  L.r) 
1 T 3'L 


sinh  (P(D)| 

w~ 


D Vg*  j (r) 


(31) 


Equations  (12)  and  (30)  represent  coupled  equations  for 
interacting  surface  and  internal  waves.  In  the  Appendix  we 
show  how  to  construct  a system  of  Hamilton's  equations  equiva- 
lent to  these.  We  shall  just  quote  these  equations  here,  using 

action-angle  variables. 

To  do  this,  we  write  for  the  surface  waves 


4>  (r , t)  = £ 
s ~ k 


2g  J, 


p w, 

. o k 


i/2 


cos(k-r  - 0.) 


C (r , t)  = ~ L 
k 


2“k  Jk 


L Pog 


1/2 


sin(k*r  - 0^) 


(32) 


Here  w = (gk) 1/2  is  the  angular  frequency  of  a surface  gravity 
}c 

wave  of  wavenumber  k in  deep  water.  The  quantity  Jk  is  the 
action  per  unit  area  of  ocean  surface  for  mode  k and  ©k  is  the 
corresponding  canonical  angle  variable. 

The  corresponding  expansions  for  the  internal  waves 


a 


2fl(j,L)  J(d,L) 


r 

cos 

L-r 

U ~ 

Here , 


2n(j,L)  J ( j , L) 


pQg  (D/L) 


T ( j ,L)  = [g'  (j,L)L2  sinh2(LD^2(j,L)] 


j(j,L)  and  0(j,L)  are  the  corresponding  action-angle  variables 
for  internal  waves  and  g'(j,L)  is  given  by  the  integral  equation 
(28)  . 

The  Hamiltonian,  referred  to  a unit  area  of  ocean,  is 
as  obtained  from  Eq.  (A. 13)  in  the  Appendix, 


H = H + H_  + V 
s I 


-l 


Jk  uk 


^ J{  j ,L)  f)<j,L) 

jrL 


- - L Wl  j.I.)  [jk  Jn  J(j.L)]  1/2 


x sin 


fe.  - 0 -0  (i  .d! 


The  quantity  G here  is 


G(k,n;j,L)  = 


2 /2p" 


fl1/2 (j ,L) 


T(jfL) 


(37) 


Also, in  writing  Eq.  (36)  we  have  dropped  as  unimportant  some 
raoidly  oscillating  terms  involving  sin  £©k  + 0n  + • 

Hamilton's  equations  are 


9H  * _ _9H_ 

30k  ' 0k  - 3Jk 


etc. 


(38) 


The  angle  variables  in  these  equations  occur  in  the  combination 


4»(k,j,L>  = 0k  - 0k_^  - 0(j  ,L)  (39) 

£note  that  the  <j>  used  here  is  not  of  Eq.  (4)j  and  the  resulting 
equations  of  motion  are 


cos<{>  (k , j ,L) 


’ - G(k+L,j,L)[jk+L  JR  J(j,L)j1/2  cos4>  (k+L,  j ,L)j 

}(j,L)  = - V 2 G(k,j,L)fjv  JV_T  J(j,L)l1/2 


cos4>  (k,  j ,L) 


<Mk,j,L)  = [wk  - w|k_Lj- 

~ Z \G(k,  jiM)[Jk_M  J(j',M)/Jk]1/2  sin<|>(k,j,/M) 

j'/M  ^ 

+ G(k+M,  j ; M)  [Jk+M  J{j# /M)/Jk  sin<J>  (k+M,  j ' ,M) 

-G(k-L,j',M)  [jk_L_MJ(  j/,M)/Jk_L]1/2  sin<j>  (k-L,  j ' ,M) 

-G(k-L+M#j;M)[jk_L+M  J(  j ' 'W)/Jk_L]1/2  sin*  3 # ,M)j 

+ ^ G(p,j,L)[jp  Jp_L/J(j,L)]  1/2sin<j)(p,j,L)  . (40) 

P 

In  these  equations  we  have  written 

G (k , j , L)  = G(k,k-L;j,L)  . (41) 

The  practicality  of  these  equations  of  motion  is  demonstrable 
in  that  integration  of  the  coupled,  first  order,  ordinary 
differential  equations  has  been  achieved  numerically  to  high 
accuracy  with  the  simplest,  first  order  Euler  scheme. 

In  the  absence  of  nonlinear  coupling,  the  Jk  and  J(j,L)' 
are  constant  and  0k  = &kt  + initial  value,  0(j,L)  = H(j,L)t  + 
initial  value,  as  is  clear  from  Eq.  (40).  We  see  then  from 
Eqs . (32)  and  (33)  that  our  choice  of  variables  is  such  that 

each  linear  mode  corresponds  to  a progressive  wave  with 
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wavenumber  vector  k,L,  etc.  The  expansions  (32)  and  (33)  may 
easily  be  modified  to  describe  standing  waves. 

We  conclude  this  section  by  noting  that  both  H and  the 
sum  of  the  J^'s  are  constants  of  the  motion  for  our  system. 

The  total  action  being  constant  is  a Manley-Rowe  condition 
[see,  for  example,  Sturrock  (I960)].  The  important  implica- 
tions of  these  two  constants  of  the  motion  for  energy  exchange 
between  two  wave  systems  have  been  discussed,  for  example,  by 
Cohen  et  al.  (1972). 
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III.  APPLICATIONS  TO  LABORATORY  EXPERIMENTS 

In  the  experiments  of  Lewis,  Lake,  and  Ko  (1974)  an 
internal  wave  was  mechanically  driven  at  a fluid  interface 
near  the  bottom  of  a wave  tank.  A small  amplitude  surface 
wave  was  also  mechanically  generated  and  the  modulation  by 
the  internal  wave  of  the  surface  wave  height  and  slope 
recorded.  The  wavelengths  of  the  surface  waves  were  so  selected 
that  their  group  velocities  were  nearly  equal  to  the  phase 
velocity  of  the  internal  wave. 

A theoretical  calculation  using  perturbation  theory 
was  made  by  Lewis  et  al.  (1974)  . This  gave  rather  good 
agreement  with  their  observations,  but  failed  to  show  a sat- 
uration at  the  larger  modulation  amplitudes. 

To  apply  our  theory  to  this  experiment  we  use  the  thin 
thermocline  model  ^Eq.  (3)]  and  adapt  the  choice  of  amplitude 
excitations  ^Eqs.  (32)  and  (33)J  to  correspond  to  the  standing 
waves  generated  in  the  tank.  For  the  thin  thermocline  model 
only  the  lowest  internal  wave  mode,  corresponding  to  j = 1, 
occurs.  Equations  (40)  were  integrated  with  the  initial 
condition  that  both  internal  and  surface  waves  were  simple 
standing  waves  at  time  t = 0.  Surface  wave  sidebands  developed 

as  a result  of  the  internal  wave  interaction  and  these  caused 
an  envelope  modulation  of  the  surface  wave  amplitude  The 
growth  in  time  of  this  modulation  was  interpreted  as  a growth 
down  the  tank,  assuming  that  the  modulation  was  convected  at 
the  surface  wave  group  velocity.  The  resulting  modulation  as 


i 

\ 


i 


I 
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a function  of  distance  down  the  tank  is  compared  with  the 
observations  of  Lewis  et  al.  (1974)  in  Figure  (1).  The  tendency 
of  the  modulation  to  saturate  at  about  22%  should  be  noted. 

In  the  above  experiment  of  Lewis  et  al.,  the  surface 
wave  amplitudes  were  not  large  enough  to  significantly  modify 
the  internal  wave.  Another  series  of  tank  experiments  to  study 
explicitly  the  generation  of  internal  waves  by  surface  waves 


were  performed  by  Joyce  (1974). 

For  our  present  purposes,  the  experiment  labeled  E4  by 
Joyce  is  of  primary  interest.  In  this  experiment  Joyce  created 
a fluid  interface  with  a sharp  density  gradient  approximately 
30  cm  below  the  surface  of  a 90  cm  deep  tank.  Two  standing 
waves  were  mechanically  driven  on  the  tank  surface  and  allowed 
to  interact.  The  experimental  parameters  were  selected  such 
that  the  interaction  of  the  two  standing  surface  waves  would 
drive  an  internal  wave  as  the  third  member  of  a resonant  triad. 
The  frequencies  and  amplitudes  of  the  two  standing  waves  were, 
respectively,  18.2003  and  17.7225  radians/sec  and  0.377  cm 
and  0.268  cm.  The  growth  of  the  resonant  internal  wave  was 

observed. 

A theoretical  calculation  of  the  resonant  wave  interaction 
using  perturbation  theory  was  performed  by  Joyce  (1974) . This 
agreed  well  with  his  observations  at  earlier  times,  but  failed 
to  show  saturation  at  larger  internal  wave  amplitudes. 

To  apply  our  theory  to  the  Joyce  experiment,  we  can 
again  use  the  thin  thermocline  model  of  Eqs.  (3).  Viscous 
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damping  terms  must  be  added  to  the  Eqs.  (40)  for  and  J(1,L), 
as  described  by  Joyce  (1974)  . The  results  of  our  calculation 
arf  shown  in  Figure  2,  where  they  are  compared  with  Joyce's 
measurements.  The  solid  curve  in  Figure  2 was  obtained 
assuming  that  the  relative  phase  <J>(k,l,L)  was  "locked,"  having 
the  value  tt.  The  dashed  curve  results  from  integration  of  the 
complete  set  of  Eqs.  (40)  . The  tendency  for  saturation  to 
occur  is  seen  to  be  more  pronounced  for  the  exact  integration 
of  Eqs.  (40) . 

The  "locked  phase"  approximation  is  often  valid  with 
fair  accuracy  in  resonant  triad  interactions.  Reference  to 
the  third  of  Eqs.  (40)  indicates  that  this  approximation  will 
fail  catastrophically  when  the  magnitude  of  the  frequency 
mismatch  becomes  greater  than  the  nonlinear  coupling  strength. 


23 


IV.  APPLICATIONS  TO  OCEAN  WAVES 

In  this  section  we  shall  relate  the  formalism  developed 
in  Section  II  to  ocean  wave  phenomena  and  discuss  some  con- 


sequences of  Eq.  (40) . 

The  power  spectrum  H^fk)  of  ocean  surface  waves  has  been 
reviewed  by  Phillips  (1966,  pg  75).  This  spectrum  is  so  nor- 
malized that 


£(-nr-)  Vi1  m J*2 kV>5>  = • 


(42) 


Here  represents  an  ensemble  average  over  many  specific 

realizations  of  the  ocean.  For  simplicity,  we  assume  that 

is  independent  of  position  on  the  ocean  surface.  In  the 
writing  of  Eq.  (42)  we  have  replaced  a sum  over  discrete  modes 
in  the  surface  region  of  area  Aq  by  an  integral  over  wavenumber. 

The  corresponding  internal  wave  power  spectrum  H'j(j,L) 
has  been  discussed  by  Garrett  and  Munk  (1972a,  1975) . The  indi- 
cated arguments  of  are  mode  number  j and  wavenumber  L.  This 
spectrum  is  so  normalized  that 


t A2l  -t  <s->  - <=2> 

j=l  J j=l 


(43) 


which  is  here  considered  to  be  independent  of  r = (x,y)  and 
to  be  evaluated  at  z = -D,  the  thermocline  depth. 
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These  power  spectra  may  be  related  to  the  action  variables 


r 


I 


■ 


using  Eqs  (32)  and  (33) . We  first  observe  that 


<v)>  = y_ [ (gpo)_1  J 


k “k 


= Ys  [Pog,(j'L)]  JU»L)  fi(j,L) 


Therefore, 


Jk  = («co/“k) 


(27T) 


^S  (~) 


J(j,L) 


f g'(j,L)p  /fl(j,L)l 

(2tt)  2 
A 

L ° J 

o 

1=3  — 

VjtjfL)  . (44) 


The  equilibrium  spectrum  suggested  by  Phillip's  (1966) 
is  of  the  form. 


(k)  = 4xl0-3  k'4  Ge<k,B)h(k)  , k > k > k 

S - S tj  vj 


h (k)  =1  for  ks  > k > kQ 


= 0 for  k 

A 

*■ 

0 

or 

; (1966)  chose 

Gs(k.B)  = | 

for 

A.  A 

k*w 

A A 

= 0 

for 

k*w 

(45) 
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(46) 


where  w is  a unit  vector  in  the  direction  the  wind  is  blowing. 
The  spectrum  suggested  by  Garrett  and  Munk  (1975)  is  for 
internal  waves  with  mode  corresponding  to  j=l: 

( j , L)  = 8 . 0X10~3  l/  [(l2  + 1.14*10"9j2)  ( j + 6)  2 * 5]  . (4 

Were  the  nonlinear  interactions  responsible  for  energy 
exchange  among  the  linear  modes  sufficiently  strong,  we  might 
expect  on  statistical  grounds  equipart ion  of  energy  with 

o>kJk  = ft(j,L)J(j,L)  = const,  all  k and  L 

This  condition  is,  of  course,  grossly  violated  because  dissi- 
pation and  generation  mechanisms  overwhelm  the  tendency  to 
equipartition  the  energy.  The  energy  per  mode  per  unit  area, 
a>kJk  and  (1  ,L)  J ( 1,  L)  , is  shown  in  Fig.  (3)  as  a function  of 
wavenumber,  using  Eqs.  (44),  (45)  and  (47),  and  assuming  a 
"rectangular  ocean"  of  dimensions  100  km  * 100  km  and  the 

-3 

Garrett-Munk  ocean  model  (2)  with  D = 100  m and  6p/pQ  = 10 
In  the  remainder  of  this  paper  we  shall  study  the 
implications  of  Eqs.  (44)  for  the  generation  of  internal  waves 
by  surface  waves.  We  can  assume  that  appreciable  exchange  of 
energy  will  occur  only  among  "resonant  triads"  satisfying  the 
condition  that  Jsee  the  third  of  Eqs.  (40)  .j 

Aw  = fu).  - a)  h _i  - Q(i,L)l  is  small.  (48) 


t 

L I 


I 


Wl 


I f 

I 

L i 


In  practice  this  condition  is  that  Aw  be  not  significantly 
larger  than  the  other  term  on  the  right  in  the  third  equation 
of  Eq.  (40) . 

We  first  describe  numerical  integration  of  Eqs.  (40) 


for  three  simple  examples.  For  each  of  these  we  have  taken 


v-3 


D = 100  m and  6p/pQ  - 1°  and  use  the  thin  thermocline  model 


of  Eq.  (3) . A single  internal  wave  mode  of  wavenumber  L 
interacts  with  three  surface  wave  modes,  i,  2,  3,  of  respective 
wavenumbers 


ki  , 


h = 5i  - i 


k3  = k2  - 2L 


(49) 


The  initial  phases  <j>  , at  t = 0,  were  all  set  equal  to  zero 
in  Eqs.  (40) . 

For  the  first  example  we  have  suppressed  mode  3 in  Eq. 
(49)  and  have  taken  (in  polar  coordinates)  the  values  listed 
under  I in  Table  I corresponding  to  exact  resonance.  The  rms 
amplitudes  of  the  three  modes  (i.e.,  ) are 

shown  in  Fig.  (4)  as  a function  of  time.  The  initial  decrease 
in  amplitude  of  the  internal  wave  mode  is  a consequence  of 
our  choice  of  initial  condition  that  the  phases  <J>  vanish  at 
t = 0.  As  J(1,L)  reaches  its  minimum  value,  the  phase  <j>(k,l,L) 
switches  to  the  value  tt,  and  J(1,L)  begins  to  increase. 


Although  not  shown  in  the  figure,  the  internal  wave  amplitude 

>4 


reaches  a maximum  of  10.9  m at  1.7x10  sec  and  then  begins  to 
decrease. 
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difference  between  the  surface  inodes  taken  two  at  a time  are  listed. 


I 


For  the  next  example,  we  have  chosen  the  modes  corres- 
ponding to  column  II  in  Table  I.  The  frequency  differences 

are  also  listed  and  we  observe  that  Aw  (=  w - w,  ) = ft(l,L) 

12  kl  k2 

The  resulting  rms  wave  amplitudes  are  shown  as  a function  of 

time  in  Fig.  (5).  The  amplitude  of  mode  3 is  seen  to  be 

nearly  constant  (to  within  less  than  one  percent) . The  reason 

for  this  is  that  it  is  too  far  from  resonance  (but  not  far 

enough  to  completely  lose  resonant  effects) . 

tor  our  third  example  we  have  chosen  the  modes  with 

wavenumbers  listed  in  column  III  of  Table  I with  the  frequency 

differences  as  shown.  Note  that  again,  Aw^  = (2(1, L).  The 

resulting  rms  wave  amplitudes  are  shown  as  functions  of  time 

in  Fig.  (6).  The  frequency  difference  Aw„_,(  = w - w ) 

23  k2  . k3 

is  now  close  enough  to  resonance  that  significant  coupling  to 
mode  3 occurs  after  the  mode  2 amplitude  has  reached  about  4 m. 

Our  third  example  suggests  the  following  generalization. 
If  we  write 

~n  = ~1  “ (n-l)L,  n = 2,  3,  ...  , (50) 

and  assume 

nL  <<  k^,  all  n's  considered, 


then 
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“k  - “k  ' ?k.'i  ■ 

n-1  n ~1 


where  C,  is  the  group  velocity  of  mode  1.  Now  if 
~1 


ft(l,L)  = C.  -L 
-1  ~ 


we  satisfy  the  resonance  condition  between  the  internal  wave 
and  adjacent  pairs  of  surface  waves.  In  this  case,  a "cascade" 
process  of  energy  exchange  can  result,  just  as  has  been  noted 
for  laser-plasma  wave  couplings  [ See  Cohen,  Kaufman,  and  Watson 
(1972)]  . 

We  next  study  the  transfer  of  excitation  from  a sharply 

collimated  spectrum  of  surface  waves  to  internal  waves.  By 

"sharply  collimated,"  we  mean  that  the  function  6 (k,$)  in 

Eq.  (45)  vanishes  unless  3 = 0 > a constant.  The  excitation 

o 

of  a single  internal  wave  mode  (j,L),  having  a very  small 

action  variable  J(j,L)  will  be  considered.  The  energy  source 

for  driving  this  internal  wave  mode  will  be  a set  of  large 

amplitude  surface  waves  having  wavenumbers  p^,  • •••  an<* 

corresponding  action  variables  J°  , J°  , ....  The  set  of 

'Pi  -P2 

waves  consisting  of  the  third  members  of  resonant  triads  are 

surface  waves  of  small  amplitude  having  wavenumbers  - L 

and  respective  action  variables  (i=l,2,...).  We  shall 

-i 

assume  that  the  large  amplitude  quantities  J°  are  constant 

?i 
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during  the  time  of  interaction  [The  approximation  of  parameteric 
amplification,  see  for  example,  Nishikawa  (1968) J.  The  "locked 
phase"  approximation  will  also  be  used,  so  Eqs.  (40)  now  read 


J*.  = 2 G(p.,j,L)  |j°  J(j,L)J 


1/2 


J(j,L)  = 


= V 2 G(p.,j,L)  [ J°  Jk  J ( j , L) 

i " 1 ~i  -i  * 


1/2 


cost})  (p^ , j ,L)  = -1,  all  modes. 


(52) 


We  look  for  a solution  to  these  equations  of  the  form 


Jk  = C.^  exp  (at) 


J ( j , L)  = D exp  (at) 


(53) 


and  find  that 


a2  = 


Z [2  Jp.  • 


Using  the  first  of  Eqs.  (44)  and  transforming  from  the  discrete 

o 

to  the  continuous  system  we  can  re-express  a as  an  integral 
over  the  surface  wave  spectrum: 


a2  =J  d2k  ^(4  gpo)/wkJ  p(k,j,L)J  2 ^(k) 


(54) 
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use  the  assumed  condition  that 


To  simplify  Eq.  (54)  we 
L <<  k and  Eqs.  (37)  and  (41)  to  '.’rite 


A o 2Q 

G lk,  j ,L)  = - 

p 

o 


,L)  [k-L  - i L21 

T2(j,L)  I ' 2 J 


Now  the  resonance  condition 


f = w I k-L  I ' f E Wk  ‘ Q(j,L)  ' 


(55) 


lets  us  write 


k-L 


-H 


k2  + L2 


f4/g2 


(56) 


or 


[k-L  - |L2j2  = 4k3  fi2  ( j ,L)  /g 


Collecting  these  results,  and  using  the  spectral  representa- 
tion (45) , we  can  write  Eq.  (54)  in  the  form 


a 


2 


efl5(i,L) (4xl0~3) 
g'(j,L)  L2  sinh(LD)_ 


dB  Gc(k,B) 

o 


; 


(57) 


i 


The  observations  of  Tyler,  et  al.  (1974) , suggest  that 
the  spectral  angular  distribution  Gg(k,B)  [Eq.  (45)]  is 
strongly  peaked  near  the  long  wavelength  cut-off  and  in  the 
direction  the  wind  is  blowing.  If  a modification  of  the 


Phillips  spectrum  (45)  is  desired,  (for  example,  to  describe 
a swell)  an  appropriate  function  h(k)  may  be  used  in  Eq.  (57) 


Now,  Eq.  (55)  implies  that  for  k = kQ;  k and  L are 
nearly  perpendicular , so  we  might  anticipate  that  a well  colli- 
mated beam  of  surface  waves  would  generate  internal  waves 
propagating  at  nearly  right  angles  to  the  direction  of  the 
surface  waves.  Apel  et  al.  (1975a)  have,  indeed,  reported 
what  seems  to  be  internal  waves  generated  by  a large  swell 
and  propagation  at  nearly  right  angles  to  the  swell. 

To  model  this,  we  take  Gs  ( k , B ) = 6(B)  and  evaluate 
the  integral  in  Eq.  (57)  as 

d6  Gs (k , 3 ) = , (58) 

c 

where  Vc  is  a "characteristic  velocity." 

+ 

We  have  taken  V = 20  m/sec  and  evaluated  a for  the 
first  four  modes  (j  = 1 , . . . 4 ) using  the  Garrett-Munk  ocean 
model  Eq.  (2).*  The  results  for  1/a  in  seconds  are  shown  in  Fig. 
(7)  as  a function  of  internal  wavenumber  L for  a mixed  layer 
depth  of  100  m.  The  solid  curves  correspond  to  6p/pQ  = 0,  the 
dashed  curves  to  6p/pQ  = 10-3  at  the  thermocline. 

The  lowest  order  internal  wavemodes  are  excited  much 
more  strongly  than  are  the  higher  modes.  Since  the  lower  modes 
are  concentrated  nearer  the  ocean  surface  this  is  not  unexpected. 

tSince  a varies  as  l//v^  the  curves  in  Fig.  (7)  are 
easily  re -scaled  for  other  values  of  V . 

*The  WKB  approximation  was  used  to  find  the  eigenmodes  (15) 


Some  insight  into  the  validity  of  the  "locked  phase" 
and  "parametric  amplification"  approximation  will  be  found  in 
the  next  section,  where  numerical  integration  of  Eqs.  (40)  is 
studied. 
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V.  GENERATION  OF  INTERNAL  WAVES  FROM  A SATURATED  SURFACE 
SPECTRUM 

In  this  section  we  shall  study  the  generation  of  a 
specific  internal  wave  mode  (j,L)  by  a saturated  surface  wave 
spectrum.  For  this  we  use  the  expression  (45)  with  G (k,0) 
given  by  Eq.  (46) . Since  the  spectrum  (45)  is  isotropic  over 
a hemisphere,  the  generated  internal  wave  spectrum  will  be 
approximately  isotropic.  A set  of  interacting  triads  will 
have  surface  wavenumbers  p^  and  = p^  - L,  i = 1,  2,  .... 

To  simplify  notation,  we  write 


ft  = ft ( j , L)  , = G(p^,j,L)  , 

Q*  S (()  (p  j ,L)  , J = J ( j ,L) 

X ^ ~ ^ 


Then,  for  the  set  of  interacting  triads  Eqs.  (40)  can  be 
written  in  the  form 
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1 
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n 


sin<|i 


n 


(59) 


n 


The  resonance  condition  (55)  can  be  satisfied  for  d.  < p , 

ri  c 

the  wavenumber  at  which  the  surface  wave  group  velocity  is 
equal  to  the  internal  wave  phase  velocity: 


pc  = (gL2)/(4ft2) 


(60) 


Thus,  we  include  only  those  modes  in  Eqs.  (59)  for  which 
pc  > Pi  > kQ  [e<3-  (45)].* 

Exact  resonance  is  not  required  in  Eqs.  (59)  for  sig- 
nificant transfer  of  excitation.  To  take  account  of  this,  we 
set 


wp  - u>k  - ft  = - 6f  , 


where  6f  represents  the  frequency  mismatch.  Now,  Eq.  (56) 
becomes,  for  a given  wavenumber  p, 

p-L  = pL  cosB  = j[k?  + L2  - (f  + 6 f ) 4/g2  ] . 


At  exact  resonance,  corresponding  to  <$f  = 0,  the  angle  6 is 


*For  our  numerical  examples  presented  later  it  can  be 
verified  that  pc  lies  below  the  range  of  capillary  wavenumbers, 
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sin8r  = 1 - p/p. 


For  6f  / 0,  we  have  8 = 6r  + <5  3»  where 


66  = (2  w 6f)/(gL  sin6r) 


The  mode  index  "i"  in  Eqs.  (59)  thus  labels  a magnitude  and 
a direction  for  each  p^. 

Integration  of  Eqs.  (59)  is  simplified  by  introducing 
dimensionless  variables,  as  follows: 


J = J Z. 
?i  ?i  1 


= Xi  ' 


JG  v , 


where  £see  Eqs.  (44) j 


J°  = ) f^l  '•’=<!?;>  = 


. 


Here  Hq  has  the  dimensions  of  length  and  is  to  be  determined 

by  the  condition  (65)  below.  A dimensionless  time  x is  intro- 
duced by  the  relation 


= a-1x  , 


where  a is  defined,  as  in  Section  IV,  by  the  equation 


“2  = I [2  Gi]  2 j°  • 


(64) 


The  constant  Hq  is  then  defined  by  the  relation 


1 

1 


2 Gi  Jo  = a(pi/Pc) 


(65) 


This  takes  proper  account  of  the  dependence  of  G.  on  p. . 

The  scaling  (63)  for  the  action  variables  implies  that 
for  the  assumed  saturated  surface  wave  spectrum  the  initial 


values  of  the  Xi  and  Z i should  be  near  unity. 


To  evaluate  a , we  can  start  with  the  expression  (57) 


If  we  make  the  arbitrary,  but  reasonable,  assumption  that 


6f|  < a , 


then  Eq.  (62)  can  be  used  to  give 

P,B)  = (4Wp0t)  / (irgL  sin&r) 


■1 


1 


1 
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r 


r 


Integration  over  wavenumber,  with  the  upper  limit  pc. 


(the  lower  limit  k can  be  taken  to  be  zero  for  the  present 


integration)  gives 


a = 


[4  *10~31 


^9*  ( j ,L)  L sinh (LD) 


(66) 


The  adopted  scaling  of  variables  in  Eqs.  (59)  leads  to 
coefficients  of  order  unity,  removing  the  coupling  strengths 


G.  and  the  dimensional  quantities  J°  and  J . Writinq 
j-  P^  o 


~ X*  ~ = 6px  6pv  = p6p'60 


lets  us  carry  out  the  mode  number  sums  as  a numerical  integration. 

Equation  (59)  was  numerically  integrated  with  the  Z.^ 
and  Xi  initially  equal  to  unity  (changing  initial  conditions 
showed  this  to  be  uncritical) . Nine  values  for  8 and  eight 
for  p^  were  used  making  a total  of  seventy -two  modes  in 
Eq.  (59)  . ^Integration  with  fewer  modes  seemed  to  indicate 
little  sensitivity  to  the  number  of  modes  used.J 


A characteristic  growth  time  tg.  for  the  internal  wave 


mode  was  defined  as 


■ </M 


evaluated  at  that  time  t that  maximum  power  was  being  delivered 


i 


1 


45 


to  the  internal  wave  [i.e.,  Y was  a maximum j .*  It  was  found 
that 


= 0.3l 

this  corresponds  to  a dimensional  time  t = t , then; 

g 

tg  = 0.31  a 1 . (1,7) 

It  might  be  noted  that  negligible  transfer  of  excitation 
occurred  for  6f  > 10  a.  |The  arbitray  choice  1 6 f | 5 a used  in 
evaluating  a is  of  course  compensated  in  the  numerical  value 
f°r  Tgi  since  is  independent  of  this  choice  .J 

Our  dimensionless  scaling  of  Eqs . (59)  was  of  course 

motivated  by  the  simple  model  calculation  of  the  last  section. 
It  is  interesting  to  note  that  our  elaborate  numerical  calcu- 
lation was  very  close  to  the  predictions  of  this  model  based  on 
the  locked  phase  and  parametric  amplifier  approximations. 

To  evaluate  tg , using  Eqs.  (66)  and  (67)  the  Garrett- 
Munk  ocean  model  ^Eq.  (2)j  and  the  WKB  approximation  were  used. 
The  resulting  values  of  the  growth  time  are  shown  in  Figs.  (8) 

- (11)  for  a range  of  parameters. 

In  Fig.  (8)  we  show  t (in  seconds)  as  a function  of 

y 

*For  Y well  below  its  saturation  value  the  choice  cf 
time  at  which  to  evaluate  t was  not  very  critical. 
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the  mixed  layer  thickness  D for  the  modes  j = 1,  L = 0.015 
and  0.005  m 1 and  with  &P/PQ  = 0* 

In  fig.  (9)  we  show  t as  a function  of  6p/pQ  for 
D = 100  m,  L = 0.01  m 1,  and  j = 1,2,3. 

Finally,  in  Figs.  (10)  and  (11)  we  show  t as  a function 

-3 

of  L for  j = 1,2,3,  D = 100  m,  and  Sp/pQ  = 0,  10  . 

The  weak  excitation  of  the  high  order  modes  is  a con- 
sequence of  their  extending  to  greater  depths.  This  greater 
depth  increases  the  value  of  g'(j,L)  ^Eq.  (28)J and  thus 
reduces  the  coupling  strength  ^This  is  seen,  for  example,  in 
Eq.  (66) j . The  very  weak  coupling  between  surface  and  internal 
waves  found  by  Kenyon  (1968)  is  related  to  this  observation. 
Kenyon  used  a constant  thermocline  model,  for  which  the 
internal  wave  extends  to  the  bottom,  leading  to  a very  large 
g'(j,L)  for  all  modes.  Since  the  coupling  interaction  occurs 
near  the  ocean  surface,  the  deeper  the  internal  wave  extends, 
the  less  effective  coupling  appears  to  be.  Our  conclusions 
are  consistent  with  those  of  Joyce  (1974)  , who  noted  that  his 
two  layer  model  gave  much  stronger  coupling  than  that  found 
by  Kenyon. 
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Vi.  VISIBILITY  OF  INTERNAL  WAVES  ON  OCEAN  SURFACE 

In  the  preceeding  sections  the  growth  of  an  internal 
wave  was  examined  under  a variety  of  oceanographic  conditions 
and  the  growth  rate  calculated  as  a function  of  the  environmental 
parameters.  Our  attention  was  focused  on  the  energy  transfer 
from  surface  to  internal  waves  via  the  resonant  triad  interaction 
and  the  modification  of  the  surface  wave  spectrum  was  not 
discussed.  In  this  section  we  examine  this  modulation  of  the 
surface  wave  spectrum  by  the  internal  wave  and  discuss  the 
"visibility"  of  the  internal  waves  from  such  surface  modifications. 

The  modulation  of  the  surface  wave  spectrum  induced 
by  travelling  internal  waves  is  visible  to  both  the  eye  and 
radar  (Ziegenbein,  1970;  Perry  and  Schimke,  1965?  Gargett  and 
Hughes,  1972).  Lafond  (1962)  suggested  that  the  internal  wave 
redistributes  extant  organic  material  on  the  ocean  surface 
leading  to  variable  slope  and  reflectivity  properties  of  the 
surface  which  are  phase  correlated  with  the  internal  wave. 

Gargett  and  Hughes  (1972)  brought  this  interpretation  into 
question  by  noting  the  directional  preference  of  the  surface 
wave  modifications  visible  in  the  Straits  of  Georgia,  i.e., 
the  surface  gravity  waves  have  long  crests  parallel  to  the 
internal  waves.  Such  a phenomenon  is  inconsistent  with  a 
surface-film  mechanism  for  producing  "slicks,"  but  not  with  a 
direct  mechanical  interaction  such  as  surface  wave  scattering 
from  the  surface  currents  produced  by  the  internal  wave. 
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The  analyses  of  the  mechanical  interaction  between 
surface  and  internal  waves  Ball  (1961);  Gargett  and  Hughes 
(1972);  Nesterov  (1972);  Lewis  et  al.  (1974)  and  Thomson 
and  West  (1975)  assume  the  internal  wave  to  be  of  sufficiently 
large  amplitude  that  the  energy  transfer  is  from  the  inter- 
nal to  the  surface  waves,  i.e.,  no  internal  wave  generation. 

To  examine  the  modulation  of  the  surface  wave  spectrum 
induced  by  an  evolving  internal  wave  spectrum,  using  a numerical 
integration  of  tne  Hamiltonian  equations,  we  consider  the  two 
model  systems  shown  in  Table  II.  In  the  first  system  we 
include  the  four  surface  modes  k^  (i=l,  2,  3,  4)  with 
amplitudes  determined  from  Eq.  (45).  The  four  internal  modes 
L.  (i  = 1,  2,  3,  4)  have  zero  initial  amplitude. 

In  Figure  12  the  rms  amplitude  of  the  internal  wave 
modes  are  depicted.  Only  mode  and  L4  have  the  appropriate 
wave  vectors  to  grow  significantly  in  this  time  interval 
(approximately  5.5  hours)  by  the  resonant  triad  interaction 
with  the  surface  wave  spectrum  selected.  The  linear  approx- 
imation to  the  internal  wave  dispersion  relation  is  tested 
by  doing  the  calculation  twice.  There  is  no  appreciable  effect 
in  the  lowest  mode  L. , but  the  second  harmonic  changes 
significantly.  The  bumps  in  the  growth  curve  are  due  to  fre- 
quency mismatch  induced  by  the  nonlinear  coupling. 

The  two  dimensional  view  of  the  ocean  surface  looking 
down  from  above  is  shown  in  Figure  13  where  contours  of  equal 
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In  column  I the  wavevectors  for  four  surface  modes 
...,k4  and  four  internal  modes  L^,  ...,L4  are  listed  in  polar 
coordinates.  In  column  II  the  wave  vectors  for  twenty  four 
surface  modes  are  listed,  where  for  each  wavenumber  (k)  there 
are  six  angles  given  by  0^  = -.38°n,  n = 0,  1,  Also 

in  column  II  are  listed  ten  internal  modes  in  polar  coordinates 
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surface  height  are  shown.  From  column  I of  Table  II  we  see 
that  the  surface  waves  are  .-.ravelling  in  a narrow  cone  from 
left  to  right  and  from  Figure  13  they  produce  a weak  inter- 
ference pattern  in  approximately  the  y direction  at  time  t = 0. 
In  Figure  14  the  negligible  initial  amplitude  of  the  internal 
wave  field  is  shown  as  contours  of  equal  height,  along  with 
the  direction  of  the  modes  included  in  the  calculation. 

Figure  15  depicts  contours  of  equal  surface  slope  (in 
magnitude)  at  time  t = 0,  ie.,  |v?(x,t)| . The  intensity  of 
incoherent  light  backscattered  from  the  ocean  surface  is  de- 
pendent on  the  mean  square  surface  slope.  By  averaging  the 
slope  contours  in  Figure  15  over  a length  large  compared  to  the 
surface  wavelengths,  and  small  compared  to  the  internal  wave- 
length, one  obtains  a measure  of  the  'visibility'  of  the  sur- 
face modulation  pattern. 

4 

After  2x10  sec  the  internal  wave  has  developed  to  an 
amplitude  of  1.44  m.  The  contours  of  equal  height  are  a 
series  of  ridges  as  shown  in  Figure  16  In  Figure  17  the 
correlation  of  the  surface  wave  modulation  pattern  with  the 
internal  wave  amplitude  is  clearly  seen.  The  small  scale 
structure  in  Figures  (15)  and  (17)  is  the  "mottling"  discussed 
ny  Apel  et  al . (1975a)  and  arises  from  the  self  interference 
of  the  surface  waves.  This  effect  is  enhanced  by  the  growth 
of  the  internal  wave. 

An  observation  of  this  long  wavelength  modification  in 
the  surface  wave  structure  (approximately  .6  km)  was  reported 
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Figure  16 


by  Apel  et  al.  (1975b)  in  discussing  pictures  of  the  ocean  sur- 
face obtained  from  the  ERTS  satellite.  This  is  the  order  of 
one  small  scale  division  in  Figures  (13)  - (17).  The  internal 
waves  were  assumed  to  be  generated  by  tidal  currents  impinging 
on  the  continental  shelf  and  the  interaction  with  surface  waves 
produced  the  visible  surface  patterns. 

A final  example  uses  a total  of  24  surface  and  10  inter 
nal  modes  as  listed  in  Column  II  of  Table  II  to  represent  the 
ocean  environment.  In  this  case,  in  addition  to  the  surface 
waves  having  the  equilibrium  amplitudes  prescribed  by  Eq.  (45) , 
the  internal  waves  amplitudes  are  given  by  the  Garrett  and  Munk 
(1972a)  spectrum  Eq.  (47).  The  resulting  initial  surface  pat- 
tern (contours  of  equal  height)  is  shown  in  Figure  18.  The 

corresponding  contours  for  the  internal  wave  field  are  shown 
in  Figure  19. 

The  response  of  the  ocean  surface  uo  this  internal  wave 
spectrum  is  depicted  in  Figure  20,  where  contours  of  equal  wave 
slope  (magnitude)  are  shown.  The  surface  shown  in  Figure  20 
is  again  qualitatively  similar  to  the  "mottling"  of  the  ocean 

surface  observed  in  photographs  taken  from  the  ERTS  satellite. 
[Apel  (1975a) 1. 
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VII.  DISCUSSION  AND  CONCLUSIONS 

The  intent  of  this  paper  has  been  to  present  a theo- 
retical model  with  which  the  energy  exchange  between  ocean 
surface  waves  and  internal  waves  can  be  calculated.  The  Ham- 
iltonian model  is,  of  course,  limited  in  its  description  of  the 
energy  transfer  process  in  an  oceanographic  environment.  The 
nonlinear  interactions  among  surface  waves,  as  well  as,  among 
internal  waves  have  not  been  included  in  the  model.  Also,  the 
coupling  of  the  wind  to  the  air-sea  interface  has  been  ignored. 
On  the  long  time  scale  of  the  surface-internal  wave  interaction 
(characteristically  several  hours)  these  phenomena  can  signi- 
ficantly influence  developing  surface  and  internal  wave  spectra. 
However,  in  the  equilibrium  case,  when  the  generation  and  dis- 
sipation (or  transfer)  rates  of  surface  waves  are  in  balance, 
the  general  conclusions  based  on  the  analysis  and  calculations 
are  thought  to  be  useful. 

The  calculated  interaction  times  given  in  Figs.  (10) 
and  (11)  can  be  used  to  estimate  th?  rate  at  which  energy  is 
fed  to  the  internal  wave  system  by  surface  waves.  The  second 
of  the  relations  (44)  lets  us  write  the  power  per  unit  area 
delivered  to  the  internal  waves  of  mode  number  j as 

Pj  = Po/g'(j'L)  Vj'L)  d2h/tg 

We  have  evaluated  this  integral  using  the  Garrett-Munk  spec- 
trum (47)  and  t taken  from  Fig.  (10).  The  integration  was 
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restricted  to  the  range  L > 7xio~4m~1,  since  Fig.  (3)  (and  our 
numerical  calculations)  suggest  that  when  L > 7xl0~4m~1  the 
internal  wave  energy  per  mode  is  too  large  to  receive  substan- 
tial energy  from  the  surface  waves.  The  results  for  j = 1,  2 
were 

2 

P1  = 3 erqs/cm  /sec 
1 2 

p2  ~ J er9/cm  /sec  • 

the  principal  contribution  coming  from  the  range 
7xl0_4m"1  < L < 5xl0“3m-1 

These  transfer  rates  are  surprisingly  close  to  the 

2 

value  of  about  1 erg/cm  /sec  obtained  by  Bell  (1975)  from 
considerations  of  internal  tides.  They  are  also  not  very 
dissimilar  from  the  dissipation  rate  of  approximately  7 ergs/ 
cm2/sec  estimated  by  Garrett  and  Munk  (1972b) . 
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APPENDIX  THE  HAMILTONIAN 


In  this  Appendix  we  discuss  the  Hamiltonian  (35)  in 
more  detail  than  was  given  in  Section  II.  That  Eq.  (35)  is  the 
correct  Hamiltonian  may  be  verified  by  simply  observing  that  the 
equations  of  motion  (38)  are  equivalent  to  Eqs.  (12),  (30), 
and  (31). 

It  is  also  easily  verified  that  the  energy  associated 
with  the  surface  waves  is 


“s  " T J 


r 2 

-f  1/2  .2 

hr 

V<D  4>sj  + gr 

( A . 1 ) 


Canoncial  variables  may  be  introduced  with  the  relations 


=-71"  v£>  ■ 

i/2 

n(r)  = - (PQg)  Ur)  , (A. 2) 


so 


HS  - I fa2*  ["2  + 9 (®1/2q)2  J 


The  equations  of  motion  are 


(A. 3) 


Ur) 


6Qlf T = "9®  Q(r)  , 
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0HC 

Q(r)  = = n(r)  , (A. 4) 

which  are  equivalent  to  the  linearized  surface  wave  equations. 
Finally,  the  Hamiltonian  per  unit  area  of  Eqs.  (35)  and  (36) 
is 


H 


s 


= H /A 
s o 


Use  of  Eqs.  (32)  and  (A. 4)  gives  the  form  (36)  of  Hg. 

The  internal  wave  energy  is  a little  more  involved 
obtain.  We  can  calculate  the  gravitational  potential  energy 
as  follows.  Let  the  displacement  of  a fluid  element  be  given 
by  two  steps.  The  first  is  a vertical  displacement  S(z,r,t). 

The  second  is  a horizontal  displacement  to  give  the  actual 
flow.  The  second  step  does  not  change  the  gravitational  energy, 
so  we  need  study  only  the  first  (virtual)  displacement  S. 

The  density  following  the  displacement  S is 

p(x)  = p(z-S)  (l  - . (A. 5) 


The  potential  energy  is  then 


Using  Eq.  (A. 5)  and  expanding  to  second  order  in  S,  we  finally 
obtain 
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2 S2  dz 


Now,  using  Eq.  (20),  we  ha  'e 
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W A exp(iL-r)  , 
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j / L 


or 
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where 


B . = A.  _ 

j , L j # L 


This  permits  us  to  write 


(A. 6) 
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using  Eq.  (28)  and  the  relation  S( 


-D,r, t)  = ^ 


(A.  7) 


The  internal  wave  kinetic  energy  is 


S,(r,t) 
j=l  3 ~ 


,»»■  nywy*--' 


KE 


(A. 8) 


where  q represents  the  horizontal  components  of  the  fluid 
velocity.  Use  of  the  condition 


5¥  + VS  * 0 


and  Eq.  (7)  allows  simplification  to  the  form 


KE  = 


= P~Y  ]T  /d2r[g' (j,<D)  ®2  sinh2(D®)/n2(j,®)]  $ . (r) 

j = l J 3 


(A.;) 


Finally,  canonical  variables  are  introduced  with  the  relations 


Q (r)  = Py2  fl  1 t *.(-,•) 
J.  ~ o D ' ~ 7 


1/2 


ni(£}  = (Pog,)  ^(^t) 


(A. 10) 


and 


j=l 

where  the  j dependence  of  Qj  and  n x is  implicit. 

The  resulting  Hamiltonian  equations  of  motion  are  seen  to 
describe  internal  waves  in  the  linear  approximation. 

The  interaction  V in  Eqs.  (35)  and  (36)  is: 


J*2*  [y  + (nQl)2J 


(A. 11) 


V = A V 
o 


(A. 12) 


= "Ao  LJdr  n(VsQ)‘[Po  ' nT  VsQlJ  * 

The  set  of  Eqs.  (12),  (30),  and  (31)  are  equivalent  to  the  set 
of  Hamiltonian  equations 


The  above  field  variables  may  be  Fourier  expanded  over  the 

unit  area  A as  follows: 
o 


where 


Sk  = (2/Ao)1,/2  sink-r 
Ck  = (2/A0)1/2  cosk-r 


The  transformation  from  the  Q,  H,  etc.  to  the  q 
is  a canonical  transformation,  the  Poisson  bracket  relations 
for  the  new  variables  being 


qk  ' pk  6k,k'  ' 


etc. 


We  also  have 


"s  - 5 Z (p*  + “k  qk)  ' 


” \ £ p(j,L)  + !)2(j,L)Q2(j,L)j  . 


A second  canonica.  transformation  leads  to  the  action 
angle  variables  used  in  tht  text: 


J. 

qk  = (2Ao  Jk/“k)2  cos0k  ' 


= ~ /2A_  w J, \2  sinG, 


/ 


Q(j,L)  = 1^2 Aq  J(j,L)/fl(j,L)J  2 COS0(j,L) 

P(j,L)  = - ^2Aq  ft  ( j , L)  J(j,L)  j 7 sin0(j,L)  . (A. 15) 

The  set  of  equations,  Eq.  (35)  in  the  text  immediately  follow 
from  (A. 15) . 
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FIGURE  CAPTIONS 


1*  The  modulation  of  the  slope  of  the  mechanically  generated 
surface  wave  by  an  internal  wave,  as  measured  by  Lewis  et  al. 
(1974)  is  compared  with  the  numerical  integration  of  the  Ham- 
iltonian equations. 

2.  The  measurements  of  Joyce  (1974)  of  the  resonant  internal 
wave  generation  by  surface  waves  [©]  are  compared  with  the 

locked  phase"  approximation  (-)  and  numerical  integration  of 
our  Hamiltonian  equations  ( ) . 

3.  The  energy  per  mode  for  surface  and  internal  waves  is  shown 
as  a function  of  wavenumber.  The  equilibrium  spectra  of  Phillips 

(1966)  and  Garrett  and  Munk  (1975)  with  D = 100  m and  6p/p 

-3  o 

- 10  are  used. 

4.  The  root  mean  square  amplitudes  for  resonant  triad  of  Eqs. 

(49)  are  shown  as  a function  of  time. 

5.  The  root  mean  square  amplitudes  for  three  surface  and  >ne 
internal  wave  node  [Eqs.  (50)]  are  shown  as  a function  of  time. 

6.  The  root  mean  square  amplitudes  for  three  surface  and  one 
internal  wave  mode  [Eqs.  (51)]  are  shown  as  a function  of  time. 

7.  E-folding  time  for  growth  of  an  internal  wave  of  wavenumber 
L due  to  Urge  amplitude  swell  is  shown.  Growth  times  are 
shown  for  the  four  lowest  modes  (j=l,  2,  3,  4)  and  a mixed 
layer  thickness  of  100  m.  The  solid  curves  correspond  to  the 
case  that  6p/p^  = o,  the  dashed  curves  to  6p/p  = 10~^. 

8.  The  growth  time  tg  [Eq.  (67)]  is  shown  as  a function  of 
nixed  layer  thickness  for  6p/pQ  = 0 and  internal  wavenumbers 
(L)  of  0.015  m"1  and  0.005  m"1. 

9.  The  growth  time  tg  [Eq.  (67)]  is  shown  as  a function  of 
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6p/pofor  a mixed  layer  thickness  of  100  m.  , internal  wavenumber 
0.01  m 1 and  modes  j = 1,  2,  3. 

10.  The  growth  time  tg  [e,.  (67)]  is  shown  as  a function  of 
wavenumber  for  a mixed  layer  thickness  of  100  m,  6p/p  = o, 
and  modes  j = 1,  2,  3. 

11.  The  growth  time  tg  [Eq.  (67)]  is  shown  as  a function  of 
wavenumber  for  a mixed  layer  thickness  of  100  m. , 6p/p  = lo-3 
and  modes  j = 1,  2,  3. 

12.  The  root  mean  square  internal  mode  amplitudes  are  shown 

as  functions  of  time  using  the  linear  ( ) and  exact  ( ) 

dispersion  relations.  The  internal  modes  are  L = (10_3m-1, 

83°)  and  L 4 = (2xlo“3m-1,  8*ej. 

13.  Contours  of  equal  surface  wave  height  are  shown  in  the 
two  dimensional  horizonal  plane  of  the  ocean  surface  at  time 
t = 0 for  the  four  surface  modes  in  Column  I of  Table  II. 

14.  Contours  of  equal  internal  wave  height  in  the  two  dimen- 
sional horizontal  plane  at  time  t = 0 for  the  four  internal 
modes  in  Column  I of  Table  II  as  shown. 

15.  Contours  of  equal  magnitude  of  surface  slope,  ie. , |vc|, 
in  the  two  dimensional  horizontal  plane  at  time  t = 0 for  the 
four  surface  modes  in  Column  I of  Table  II  as  shown. 

16.  Contours  of  equal  internal  wave  height  are  shown  in  the 
two  dimensional  horizontal  plane  of  the  ocean  surface  at  time 

t = 2*104  sec  for  the  four  internal  modes  in  Column  I of  Table 
II. 


78 


17.  Contours  of  equal  magnitude  surface  wave  slope  in  the  two 

4 

dimensional  horizontal  plane  of  the  ocean  surface  at  t = 2x10 
sec  for  the  four  surface  modes  in  Column  I cf  Table  II  as  shown. 

18.  Contours  of  equal  surface  wave  height  are  shown  in  the 
two  dimensional  horizontal  plane  of  the  ocean  surface  at  time 
t = 0 for  the  twenty  four  surface  modes  in  Column  II  of  Table 

II. 

19.  Contours  of  equal  internal  wave  height  are  shown  in  the 
two  dimensional  horizontal  plane  of  the  ocean  surface  at  time 
t = 0 for  the  ten  internal  modes  ii  Column  II  of  Table  II. 

20.  Contours  of  equal  magnitude  surface  wave  slope  are  shown 
in  the  two  dimensional  horizontal  plane  of  the  ocean  surface 

at  time  t = 103  sec  for  the  twenty  four  surface  modes  in  Column 
II  of  Table  II. 
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